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Effect of reorientation statistics on torque response of self propelled particles 
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We consider the dynamics of self-propelled particles subject to external torques. Two models for 
the reorientation of self-propulsion are considered, rnn-and-tnmble particles, and active Brownian 
particles. Using the standard tools of non-eqnilibrium statistical mechanics we show that the run 
and tumble particles have a more robust response to torques. This macroscopic signature of the 
underlying reorientation statistics can be used to differentiate between the two types of self propelled 
particles. Further this result might indicate that run and tumble motion is indeed the evolutionarily 
stable dynamics for bacteria. 


I. INTRODUCTION 

Self-propelled particles are inherently out of equilib¬ 
rium objects that consume energy at the scale of individ¬ 
ual entities to produce persistent self generated motion. 
Models of self propelled particles have become the proto¬ 
type theoretical system used to understand the physics 
of active materials. These models are directly applica¬ 
ble to micron sized objects such as motile bacteria and 
self-phoretic colloids m- 

There have been numerous theoretical investigations 
of the collective behavior of active particles in the re¬ 
cent years. Phenomenology uncovered includes cluster 
formation, phase separation, motility induced segrega¬ 
tion ims], anomalous behavior of mechanical properties 
such as pressure miGzi and rheology nans]- Further, it 
has been established that the presence of external forces 
or confinement [Ml |M] dramatically alters the 
observed phenomenology in these systems. One striking 
example has been rectification [M1IM1 - I5n] - in which asym¬ 
metric barriers induce directed transport of self-propelled 
particles. 

Generally, self-propelled particles travel along a body 
axis u that has an intrinsic reorientation process. Two 
model classes of self-propelled particles that have been 
used extensively are Run-and-Tumble particles (RTFs) 
and Active Brownian particles (ABPs). The only dif¬ 
ference in the dynamics of the two classes is the nature 
of their reorientation process. The dynamics of an RTF 
consists of periods of straight line motion followed by 
a sudden tumble event occurring at some mean tumble 
rate a. The tumble event completely decorrelates the 
orientation and a new orientation is selected at random. 
There are several examples of motile bacteria that follow 
run-and-tumble dynamics |31j . the most famous exam¬ 
ple being Escherichia co(j. [32]. The second class, ABPs, 
reorient their body axis through gradual rotational diffu¬ 
sion and this diffusion is Brownian with rotational diffu¬ 
sion coefficient Dr- These dynamics have been observed 
in experiments of synthetic self-phoretic colloids [IHS]. 
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For times much longer than the angular reorientation 
time, that for RTFs and t 7^ for ABPs, both 

classes of particles exhibit diffusive behavior. Equating 
the diffusivities at this stage would yield identical phe¬ 
nomenology for both models. This idea has been explored 
extensively in [T1|T31 |U1 [Ml 1331 |M| where the authors find 
equivalence between the two models at the level of a drift 
diffusion equation when a = (d — l)Dr, where d is the 
spatial dimension. 



FIG. 1. (color online) Orientation probability distribution 
functions for RTFs (red/solid) and ABPs (blue/dashed) in 
the case of both a nematic aligning (left) field and a polar 
aligning field (right). RTP distributions are sharply peaked 
about the optimal direction while the ABP distribution is the 
smooth equilibrium result. 

In this work we seek to identify signatures of the re¬ 
orientation statistics on the long time behavior of self- 
propelled particles. In particular we ask the following 
question : If the swimming parameters are equivalent, 
how does the long time behavior of ABPs and RTFs dif¬ 
fer when subject to external torques? External torques 
model gradient seeking behavior exhibited by microor¬ 
ganisms such as chemotaxis 1311 Ea, and aerotaxis [35j . 
Such phenomena need not be restricted to biological 
systems, as micron sized platinum-gold rods have been 
shown to exhibit directed movement towards regions of 
higher hydrogen-peroxide concentrations |3hj . By using 
the framework of non-equilibrium statistical mechanics 
we show that run-and-tumble particles align with exter¬ 
nal torques (or gradients) much more robustly than ac¬ 
tive Brownian particles, indicating that the reorientation 
behavior of bacteria is ideal for successful gradient seek¬ 
ing. 
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II. THEORY 




and the nematic order parameter tensor Q 


Let us consider non-interacting self propelled particles 
moving in unconstrained 2D space. The statistical me¬ 
chanics of this system is given in terms of Uj t), the 
probability of finding a particle at a point r at a time 
t oriented in a direction u = cos Ox + suiOy in space. 
When the particle reorients its direction through rota¬ 
tional diffusion, the dynamics of this probability distri¬ 
bution obeys a diffusion equation of the form 

dtipir, u, t) = -V • [vuijj{r, u, t)] -f V • [DtV-ipir, u, t)] 

+ de[Drdetljir,u,t)] 

( 1 ) 

where v is the self propulsion speed of the particle, Dt is 
a translational diffusion constant and Dr is a rotational 
diffusion constant. When the particle undergoes run and 
tumble dynamics, the probability distribution obeys a 
master equation of the form 

dt'ipir, u, t) = -V • u, t)]+W ■ [DtWip{r, u, t)] 

- + -^ j duV(r,u',t) 

( 2 ) 


where a is the tumble frequency. This description is valid 
when the duration of a tumble is short compared to the 
time scales set by the self propulsion and the tumble fre¬ 
quency of the particle. We are interested in the dynamics 
of this system in the presence of uniform external torques. 
So, in the following, we consider the homogeneous limit 
of Eqs. (T][2) in the presence of a 0 dependent potential. 
The dynamical equations of interest then are of the form 


dtil}{9,t) = dg[Drdgil}{9,t)] + ^dg[dgV{9)'tp{0,t)]] (3) 

for active Brownian particles and 


dt-ipiO, t) =- ^) + -^ J dO'-ipiO', t) 


-dg[dgV{9)^{9M 


(4) 


for run and tumble particles. In the above is a rota¬ 
tional friction constant characteristic of the medium in 
which the particles move. 

We will now consider different choices for the external 


potential and construct series solutions to Eqs. (3 4) of 
the form ^(0) = oq -f cos{n9) + sin(n^ for 

different choices of external potential V {9). Other use¬ 
ful measures to characterize the response of the system 
are the first two moments of the orientation distribution, 
namely, the polarization P defined as 


Pa = J dOUatpi^O) 


QajH 



]^5ap)lp{9) 


These will be calculated as well. 


III. NEMATIC TORQUE 


First we consider an external potential of the form 
V{9) = —7COs(20). Such a field will induce a nematic 
aligning torque on the particles because of the two min¬ 
ima located at 0 and tt. In this case a good characteriza¬ 
tion of the degree of order will be given by the component 
of the nematic order parameter tensor describing align¬ 
ment in the x direction, namely Qxx- 



FIG. 2. (color online) RTFs exhibit higher nematic order than 
ABPs as shown by the nematic order parameter being higher 
for RTFs (red) than for ABFs (blue/dashed) 


In the following, the parameter n is the ratio of the 
field strength of the external potential to the friction co¬ 
efficient in units of the angular reorientation time and 
is given by k = for active Brownian particles and 

K = for run and tumble particles. Computing Qxx 

from one can show (see Appendix]^ for details) that 

for active Brownian particles, we obtain 


~ 2 Jo(«) 


(5) 


which is the equilibrium result for a thermal nematogen 
at temperature ^'"Dr in the presence of an external field. 
Here the A’s are modified Bessel functions of the second 
kind m- For run and tumble particles, the exact value of 
the order parameter turns out to be a continued fraction 
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of the form 


Q 


XX — 



( 6 ) 


As seen in Figj^ the nematic order parameter is larger 
for run-and-tumble particles. 

We can also exactly compute the orientation distribu¬ 
tion functions in the presence of this external torque. For 
active Brownian particles the distribution is of the form 
(Appendix 


m 


g/t cos(20) 

27rJo(K) 


(7) 


which is precisely the equilibrium result and is usually 
called a Von Mises distribution [551 I3S]) the analog of 
the Gaussian distribution on a unit sphere. For run and 
tumble particles we obtain the following formally exact 
series solution 


OO 

^{0) = Oo -f E a 2 n(«^) cos(2n0) (8) 

n—1 

where the coefficients are given by 
2 

02 = -Qxx and a 2 „ = a2„_4 - : n > 2 

TT ^ ' 

The orientation distribution is much more peaked around 
the minima of the external potential for run and tumble 
particles than for active Brownian particles (see Figj^, 
another manifestation of the stronger response that RTFs 
exhibit to external torques. This more robust response to 
external torques is a consequence of the run-and-tumble 
search strategy being able to sample the optimal direc¬ 
tions more efficiently than its Brownian counterpart, as 
the tumble directions are entirely random and unbiased 
by the initial orientation. 


IV. POLAR TORQUE 

In this section we consider particles subject to an ex¬ 
ternal potential of the form V{9) = — 7 C 0 s(d), i.e., a 
torque that tends to align their direction of motion along 
one direction in space. In this case, the relevant measure 
of the response of the system is the polarization. For 
active Brownian particles the polarization is found to be 
(Appendix 


Px 



(9) 


while for run and tumble particles we have 


Px = 


( 10 ) 


2 + 


1-f 


^ -f 

A ' 


where k is as defined earlier. As in the nematic case in 
section 3, the result for the active Brownian particle is the 
same as that of equilibrium. The orientation distribution 
functions for the two classes are found to be 


m 


gK COS(0) 

2ttIq{k) 


( 11 ) 


for ABPs, while for RTFs 


OO 

'0(0) = oq + y] a„(/c) cos(n0) (12) 

n—1 


where the coefficients are given by 

Qi = Fj, and Cln — 0,^—2 (n—ViK ■ ^ ^ 2 
TT ' ' 

These results are plotted as a function of the dimension¬ 
less field strength k in Fig 3a and Fig 1. Unlike in the 
nematic case, at the level of the moments of the distri¬ 
bution, the ABPs exhibit higher polar ordering than the 
RTFs even though the distribution itself is more sharply 
peaked about the optimal direction. A better measure 
of the response in this case is the percentage of particles 
that find the optimal direction. We find that in terms of 
this measure the response of the run-and-tumble particle 
is indeed much more robust (see Fig([^). 



FIG. 3. (color online) (a): The polarization as a function 
of the parameter k . The polarization for RTFs (red/solid) 
is less than for ABPs (blue/dashed) (b):The percentage of 
particles oriented within some range [tt/IS, tt/IS] as a func¬ 
tion of the parameter k . It is clear that the percentage of 
RTFs (red/square) oriented in that range is higher than ABPs 
(blue/circle). 


V. DISCUSSION 

We have shown that the response of self-propelled par¬ 
ticles to external torques is sensitive to the nature of the 
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reorientation statistics of the propulsion direction. We 
find that run and tumble particles exhibit a more robust 
response to applied torques in the following sense. When 
there is one optimal direction for the particle’s orienta¬ 
tion, such as in the polar case, the fraction of particles 
that find the optimal direction is greater for those parti¬ 
cles using a run-and-tumble search strategy, even though 
the average polarization is higher for particles undergo¬ 
ing Brownian rotational diffusion. For cases when there is 
more than one optimal direction for the particles to move 
in (such as the nematic case, where there exists 2 optimal 
directions) the run-and-tumble strategy always leads to 
better alignment as the particles are able to sample the 
different optimal directions better than through Brow¬ 
nian diffusion. These external torques model gradient 
seeking behavior exhibited by microorganisms and the 
more robust response of run-and-tumble particles might 
indicate a biological preference toward this search strat¬ 
egy- 

Even though the torque response is dramatically differ¬ 
ent based on the reorientation statistics, a number of bulk 
phenomena such as phase separation and clustering are 
identical in the two classes of particles. This dichotomy 
can be understood by noting that the two models of self 
propulsion considered here have identical correlations but 
their response is very different. While the response of 
the Brownian particle has the conventional fluctuation- 
dissipation relationship to the correlation, the run-and- 
tumble particles do not have this property (see Appendix 
C). Therefore the nature of the reorientation statistics 
becomes important when considering external perturba¬ 
tions such as confining fores or aligning torques. 

Finally we note that the exact expressions obtained in 
this work for run-and-tumble particles are given in terms 
of formal infinite series. These series solutions converge 
very slowly (see Appendix B) and any truncation at low 
orders drastically fails to capture the true nature of the 
solution. Therefore, low moment closure estimates, con¬ 
ventionally used in the literature of active particles, is un¬ 
reliable for the case of run-and-tumble particles though 
they work very well for the case of active Brownian par¬ 
ticles. 
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Appendix A: Derivation of order parameters and 
distribution functions 


In this appendix, the details of the derivation of steady 
state solutions to equations (3]|4) for different external 
torques are given. Let us begin by considering the ne¬ 
matic torque. The statistical mechanics in this case is 
given by a Fokker-Planck equation for active Brownian 
particles 

5tV^(6>,t) = de[Drdeilj{0,t)]+‘^dg[sm{29)tj}{e,t)]] (Al) 

and the following master equation for run and tumble 
particles 


dtipiO, t) =- ai>{0, t) + ^ J d9'i^{9\ t) 


|^ae[sin(26»)V’(6»,t)]] 


(A2) 


As described earlier we search for series solutions '0(0) = 
Oq + cos(n0) -I- hn sin(n0). This Fourier decom¬ 

position is equivalent to the following moment expansion 

0(0) — (0 A T • . • ) 

ZTT I 

In this form it is apparent that Qxx = f 02 - Using this 
with Eq.( AI) we arrive at the following hierarchy of equa¬ 


tions for active Brownian particles 


5*00 = 0 


(A3) 


dt^ndn,m — T SgTl(^Tn U,)(l2,|m—n| ) 

r* Drdn^n^m T ^^00*^2,m 

(A4) 

(^2,n-t-m T Sgfl{Tl m)(l2,|m—n|) 

Tl Dr9ndn,m 

(A5) 

And similarly we find the corresponding hierarchy for run 
and tumble particles. 


5*00 = 0 


(A 6 ) 


9t^ndn,m — '^^ni^2,n+m T SgTl(^Tn n.)(l2,|m—n| ) 

(A7) 

5*0n*^n,m ~ (^2,n-t-m T SgTlijl ^)f^2,|m—n|) 

(A 8 ) 

First note that all odd coefficients are zero by the require¬ 
ment the 0(0) = 0(0 + tt). By truncating at arbitrary 
0271+2 or 02 n +2 and taking the steady state solution, one 
finds that all 02 „ or 6271 couple to the zeroth coefficient 
oo or &o- Since 60 does not exist, the 6277 , vanish and by 
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using an iterative procedure, we solve for 02 and find for 
active Brownian particles 


02 — 2ao- 


(A9) 


4 + 


K 


; + ... 


while for run and tumble particles we have 

K 

02 = 2oo- 


+ 


K 


(AlO) 


i + 


g + y 


Beyond Eq. (|^ no further analytic treatment for run and 
tumble particles is available. For active Brownian par¬ 
ticles the continued fraction Eq. (A9) is of the form of a 

lUSS CO 

I ident: 

Iu{z) 
h-l{z) 


Gauss continued fraction. In general we have the follow¬ 
ing identity for the ratio of modified Bessel functions 


2u- 


2(z/ “h 1) + 


2{u -h 2) “h ' 


2(z/ -h 3) -h ■.. 

(All) 

One can find continued fraction representations for all the 
coefficients in our series solution and by using Eq.(All) 
one finds that 


0-271 — 


dri('t) 


(A12) 


thus the orientation distribution function is given by the 
exact form 

1 7 00 

= ^ + ^-rw5I^(K)cos(2n6i) (A13) 

27r 7r/o(K) 

With the help of the following identity, an example of a 
Jacobi-Anger expansion 

00 

gZcos(2e) =/g(^) + 2y]4(K)cos(2n0) 
n—1 


Eq.(Al3l is able to be summed exactly to the form of 
The same process holds for the polar aligning 
torque and repeating the calculation in that case yields 
Eq.@-@ 


Appendix B: Truncation error and perturbation 
theory 

As stated in the main text, the series solution for run 
and tumble particles is slowly converging and a low mo¬ 
ment approximation is not valid in this case. One way 


to see this is to look at the % difference in the value 
of the function when you truncate the series solution at 
successive terms (see Fig. ©)• Even in the weak field 
(k = .3), a difference of about 15% is seen when truncat¬ 
ing at the second or third terms in the series. For active 
Brownian particles, the same truncation yields a 1.4% 
difference, which shows that a low moment approxima¬ 
tion is reasonable for ABPs. To obtain a % difference of 
this magnitude for RTFs, one must truncate the series at 
n = 300. For both classes, the number of terms needed 
increases as k gets larger as should be expected. 


RTP, k=.3 



2 4 6 8 10 12 

tt of terms included in series 


FIG. 4. (color online) For RTFs a large number of terms in the 
series must be included. Above shows the percent difference 
in the value of ip{0) whether truncating at the nth term or 
nth + 1 term. Here we examine the value at 0 = 0. The first 
point on the horizontal axis represents the percent difference 
in 4’{S) by truncating the series at n = 2 and n = 3 , the next 
point represents the percent difference when truncating the 
series at n = 3 and n = 4 this is continued up to the the last 
point which represents the percent difference by truncating 
the series at n = 14 and n = 15. 


Appendix C: Relationship between Correlation and 
Response 


As another measure of the difference between the two 
classes of active particles one can compute the correla¬ 
tion and response functions. To illustrate, consider RTFs 
in the presence of a polar aligning torque. To compute 
the correlation and response functions we will need the 
average of cos(0) , i.e., the polarization given by the first 
moment of the angular orientation distribution function. 



00 

d0cos(0)[ao + y] a„(K) cos(nd)] 

n—1 


(Cl) 
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= 7rai(K) = 27rao 


1 + 


1 + 


FIG. 5. (color online) The ratio of correlation to response for 
RTFs in the presence of a polar aligning torque (Red/Dashed) 
and a nematic aligning torque (Blue/Dot-Dashed) as a func¬ 
tion of K. The constant line at 1 represents the value of this 
ratio for ABPs 


Where we used Uq — l/27r. The response function is Similary, for a nematic aligning field the relevant quan- 
giving by tity is 


R=^^(cos{e)) = \-\K^+O{n^) (C2) 

The response function measures how the order parameter 
changes in response to a changing field. The correlation 
function is given by 


C = 


(cos^{9)^ — (cos{9)^ 


(C3) 


1 

2 





02 = Oi- 

1 + 



^cos(20)^ = / d0cos(20)[ao + E <^ 2 n(^) cos(2n^)] 

(C5) 

n K 

= 7ra2(^) = 27rao - - - a- - --- = —— 

K. . K. 

— -I- - — -l- - 

i J-1 _|- 

4^1, 4^1, 

6 -e . . . g -h . . . 

Where we used oq = l/27r. The response function is 
giving by 

i? = E ^cos(26»)^ = 2 - 48^2 -b 0{k^) (C6) 

The correlation function is given by 


d-d 

2 4 


+ 0{k^) 


So up to quadratic order in n we have 


C = (cos^{29)^ — ^cos(20) 


(C7) 


C 

R 



(C4) 


1.5 

1.0 

o 

0.5 


0.0|_ 

0.0 0.1 0.2 0.3 0.4 0.5 

K 


1 

2 


-a^in) - a2(K)^; 


O 4 = 02 - 


1 + 


6 + y 


^ C = + 0{k^) 

So up to quadratic order in n we have 

^ = ^+Qk^ + 0{k^) (C8) 

i? 4 

The exact expression for both aligning torques is ex¬ 
plored numerically and shown in Fig.5. Repeating this 
process for ABPs would yield CjR = 1 irrespective of 
the value of k. 
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